Quantum order by disorder and accidental soft mode in Er 2 Ti 2 07 



M. E. Zhitomirsky, 1 M. V. Gvozdikova, 1 P. C. W. Holdsworth, 2 and R. Moessner 3 

1 Service de Physique Statistique, Magnetisme et Supraconductivite, 
UMR-E9001 CEA-INAC/UJF, 17 rue des Martyrs, 38054 Grenoble Cedex 9, France 
Laboratoire de Physique, Ecole Normale Superieure de Lyon, CNRS 69364 Lyon Cedex 07, France 
s Max- Planck- Institut fur Physik komplexer Systeme, 01187 Dresden, Germany 

(Dated: March 2, 2013) 

Motivated by recent neutron scattering experiments, we derive and study an effective "pseudo- 
dipolar" spin-1/2 model for the XY pyrochlore antiferromagnet Er2Ti207. While a bond-dependent 
in-plane exchange anisotropy removes any continuous symmetry, it does lead to a one-parameter 
'accidental' classical degeneracy. This degeneracy is lifted by quantum fluctuations in favor of the 
noncoplanar spin structure observed experimentally — a rare experimental instance of quantum order 
by disorder. A non-Goldstone low-energy mode is present in the excitation spectrum in accordance 
with inelastic neutron scattering data. Our theory also resolves the puzzle of the experimentally 
observed continuous ordering transition, absent from previous models. 

PACS numbers: 75.10.Jm, 75.30.Ds, 75.50.Ee, 



Introduction. — Exotic low-temperature properties of 
magnetic pyrochlore oxides have their origin in the un- 
derlying frustrated lattice of corner-sharing tetrahcdra, 
which is combined with diverse crystal-field effects for 
rare-earth ions [T] . Initially, the attention was focused on 
to spin-ice materials Ho 2 Ti 2 07 and Dy 2 Ti 2 07 with Ising- 
like moments [S] [3] and on Heisenberg antiferromagnets 
Gd 2 Ti 2 7 and Gd 2 Sn 2 7 with large S = 7/2 spins [IE]. 
The above pyrochlore compounds exhibit various kinds of 
ordered and unusual disordered classical magnetic phases 
naturally expected in view of large rare-earth moments. 
However, it was gradually recognized that collective be- 
havior of low-energy magnetic doublets selected by the lo- 
cal crystalline electric field (CEF) for other rare-earth py- 
rochlores should be generally described by effective non- 
Ising pseudo-spin-1/2 Hamiltonians implying significant 
quantum effects [5H5]. 

Er 2 Ti 2 07 has a special place in the family of magnetic 
pyrochlores. It is the only pyrochlore antiferromagnet 
found so far that exhibits strong easy-plane anisotropy 
with local moments confined to planes orthogonal to the 
local (lll)-axes. Erbium Titanate has been extensively 
studied in the past and is well characterized experimen- 
tally [TUl - fT?] . It orders below T c ~ 1.2 K, into an unusual 
k = non-coplanar antiferromagnetic structure shown in 
Fig. [lja) , the so called vb 2 state [TSJ, [T3] . This is rather 
surprising given that the geometry of the four local easy 
planes also allows coplanar spin arrangements such as the 
the ip3 state shown in Fig. [ljb) , with the same energy as 
the noncoplanar state. Moreover, states obtained by re- 
versing a pair of antiparallel spins in Fig. [TJb) minimize 
the dipolar energy [15] , 

Quantum order by disorder effect Q15] was suggested 
as selection mechanism for the noncoplanar spin struc- 
ture [12) . because similar thermal order by disorder se- 
lection [20] indeed occurs in a model of classical XY 
moments [T2l [2"Tl - f2"5] . However, the concomitant tran- 
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FIG. 1: (Color online) Possible q — spin structures of the 
XY pyrochlore antiferromagnet in zero applied field: ip2 state 
(a) and ip3 state (b). 



sition turns out to be first order, in clear contradiction 
with the continuous second order transition observed in 
Er 2 Ti 2 07 [TDJ[T5]. Another recent experimental result, 
which has so far remained unexplained, is the presence 
of a quasi-acoustic mode in this highly anisotropic anti- 
ferromagnet found by inelastic neutron scattering (INS) 
measurements 14J. 

In this paper we formulate a minimal spin model for 
Er 2 Ti 2 07 which accounts for this entire body of exper- 
imental results: quantum order by disorder stabilizes 
the observed ground state, with the low-energy mode in 
INS a remnant of the lifted degeneracy; while classical 
Monte Carlo simulations reproduce the continuous tran- 
sition. Using this model, we analyze theoretically var- 
ious physical effects such as the energetic selection due 
to quantum fluctuations, spin-wave results for the exci- 
tation spectrum and the nature of the finite-temperature 
phase transition studied by classical Monte Carlo simu- 
lations. We present various estimates of parameters in 
the model Hamiltonian, which has analogous form — but 
opposite sign! — to a dipolar interaction. 

Model— Free Er 3+ ions (J = 15/2, g 3 = 6/5) have 
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large magnetic moments /i = g.jfJ.B\/ J (J + 1) = 9.58/ig. 
A somewhat reduced value p, — 8.9 ps found in Er2Ti2 07 
from high-temperature susceptibility measurements [11] 
reflects a partial quenching of the angular momentum 
by CEF. According to crystal- field analysis [121 l2~4"] . 
the ground-state Kramer's doublet |±) of Er 3+ ions is 
strongly anisotropic: 

(±|J z |±)=±iA,, (-|J-|+) = (+|J+|-) = A ± . (f) 

with A ± = 6.35 and X z = 0.202. Neglecting higher CEF 
levels, one can use the operator equivalence 



J z = KS Z 



(2) 



where S a are components of an S = 1/2 pseudo-spin de- 
fined in the subspace of the lowest doublet. Note that the 
operator relations ([2| are satisfied in the local coordinate 
frame for each magnetic ion. 

In order to discuss possible forms of the effective 
Hamiltonian, it is instructive first to assume isotropic 
exchange between full magnetic moments, J, • Jj, and 
to project the nearest-neighbor (n.n.) exchange Hamil- 
tonian onto the subspace of lowest CEF-doublets with 
the help of Eq. Q [55]. Terms with Sfo are further- 
suppressed by a small factor A z /Aj_ ~ 10~ 2 and are ne- 
glected below. The remaining XY interactions are gen- 
erally characterized by two parameters [26) : 



H = 



s- 1 



,)(S 



(3) 



(ij) 



where = (r, — rj)/|rj — rj\. The first term is isotropic 
exchange between transverse components of n.n. spins 
confined to the local {111} planes. The second term rep- 
resents bond-dependent exchange anisotropy. It has a 
pseudo-dipolar form but both signs of J a are allowed 
by symmetry. Although the above projection procedure 
yields vanishing J a , the in-plane exchange anisotropy 
may be generated by virtual excitations into higher CEF 
levels [25] or into states with different occupancy on rare- 
earth ions [7] . As we shall see below, it is essential to in- 
clude J a in the effective Hamiltonian ([3]), since quantum 
fluctuations are singular in the J a — > limit. 

Ground-state selection. — The classical ground states of 
the Hamiltonian ^ with J a = satisfying the constraint 
Etet = f° rm an infinite but nonextensive set of 
states, producing planes of zero-energy modes in the ex- 
citation spectrum |22| . Finite anisotropy J a ^ almost 
eliminates this degeneracy. For J a > 0, we find that the 
classical energy ^ is minimized for a subset of k = 
states. Two of them, so called ip2 and ^3 states [15] . 
are shown in Fig. [T] while other spin configurations are 
constructed as their linear combination: 



S n (ip) = S (x„ cos (p + y n sin cp) 



1. 



(4) 



such that ip — {ir/2) corresponds to ip2 (tps), respec- 
tively. Thus, instead of an infinite number of variables, 



the classical ground-states of Eq. ^ for J a > are pa- 
rameterized by one global continuous angle, ip. 

Spin states ip2 an d ip3 form a basis of the two- 
dimensional irreducible representation E of the tetrahc- 
dral point group. They transform as ip% ~ 3z 2 — r 2 and 
-03 ~ x 2 — y 2 under cubic rotations. Consequently, the 
two spin configurations have the same mean-field energies 
for any extension of the Hamiltonian (J3J , involving only 
bilinear spin-spin interactions {e.g., further-neighbor ex- 
change or long-range dipolar interactions). The mean- 
field degeneracy may be lifted only by interactions of 
sixth order in spin components [27j . though these are 
not explicitly present in the effective Hamiltonian since 
S = 1/2. Nevertheless, such interactions may be gener- 
ated dynamically by thermal or quantum fluctuations. 

To study the effect of quantum fluctuations on the re- 
maining continuous degeneracy we use harmonic spin- 
wave theory. For this, spin operators in the effective 
Hamiltonian ([3]) are assumed to have an arbitrary spin 
S and the usual 1/S expansion is applied, substituting 
S = 1/2 in the final expressions. This standard approach 
yields the following results. 

Magnon excitation energies for an arbitrary classical 
ground-state Q can be expressed as 



.(«) 



= JS^{2 + Ja ) 2 + 2(2 +j a )X 



(n) 
k 



(5) 



where, j a = J a jJ and X^f 1 are four eigenvalues of the 
matrix Ck, 



with c± y 



t X y COS(k x ± ky)/4, 

4 



/ 





C+ 


C+ 


C xz \ 




c+ 

x y 





^xz 


C yz 




C+ 

C+ 
^xz 


C xz 





r 

^xy 

/ 


V 


c yz 


c~ 
^xy 



(6) 



and 



sin(<^ — 7r/3) sm(ip + 7r/3) — sm 2 ip , (7) 



txz 



-- sinips'm((p — tt/3) — — sin 2 (t/? + tt/3) , 
3 3 

-— sin ip sin (ip + tt/3) — ^ sin 2 (ip — tt/3) . 



The quantum correction to the ground-state energy 
normalized per spin is 



n— 1 



d 3 fc 



1 



V BZ ~ k 2 



(8) 



where Vbz is the Brillouin zone volume. AE S _ B , (ip) is 
shown in the inset of Fig.[2]for 5=1/2, exhibiting a clear 
minimum for the noncoplanar ip2 state (tp = 0,7r/3, . . .) 
for any value of j a . Hence quantum fluctuations, like 
thermal fluctuations, select the i[>2 state from the ensem- 
ble of classical ground states giving an example of quan- 
tum order by disorder. A similar calculation for isotropic 
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FIG. 2: (Color online) Inset: the ground-state energy in units 
of J for degenerate classical configurations parameterized by 
angle ip, see Eq. |4j. Main panel: the quantum spin reduc- 
tion in the noncoplanar state rp2 versus the bond-dependent 
exchange anisotropy parameter j a = J a j J ■ 



exchange j a = 0, was recently presented in [23]. However, 
as shown below, quantum fluctuations become singular 
and essentially nonperturbative in the isotropic case and 
may completely destroy the long-range magnetic order. 

Harmonic spin-wave theory also yields the leading 
quantum correction to the value of the ordered moment 
(S). Most conveniently this can be done with the help 
of a fictitious staggered field h s 0. The corresponding 
expression for SS = S — (S) is 



ss = 



-t 

8 1 



d 3 k 
Vbz 



2+ Ja + x { ; i) 



(2 + i )2 + 2(2+i )A^ 



(9) 



The spin reduction SS in the ordered ip2 state (Fig. [2| 
diverges logarithmically as j a — > 0. This is a result of 
the effective dimensional reduction to one dimensional 
behavior, due to the planes of zero-energy modes in the 
magnon dispersion in the isotropic case [55]. How- 
ever, a small anisotropy of j a ~ 0.05 already cuts off 
the singularity and restores the perturbative nature of 
the quantum correction. The measured ordered moment 
of fi = 3.2fiB in Er 2 Ti 2 07 [13] indicates a 16% renor- 
malization from the nominal value fj,g « 3.8/xb for the 
Er 3+ ground-state doublet [T2]. This gives an estimate 
of j a ~ 0.3-0.5. The large uncertainty is related to the 
unknown precision for fj,o, which cannot be measured di- 
rectly, but our analysis is consistent with a perturbation 
of considerable size. 

Excitation spectrum. — With both Hamiltonian and its 
ground state in hand, we compute the spin-wave spec- 
trum in Er 2 Ti 2 07. The ip2 magnetic structure allows six 
different domains with tp = nn/3. We choose ip — (pic- 
tured in Fig. [lja)) with spins pointing along [±1,±1,±2] 
cubic directions. Analytical expressions for magnon en- 
ergies can be derived from Eq. |5| for a few symmetry 
directions in the Brillouin zone. For k = (g,0, 0), the 



four branches are 

2 = (2+Ja) [2+JaT 2 cos |± |i Q (l+cos |)] . (10) 

Along the k = (0,0, q) line, the two dispersive branches 
are expressed as 



£k = JSy/(2+j a )j a [l±cas(q/4)] 



(11) 



while the two other modes remain dispersionless with 
£k = JS^/(2+j a )(A + j a ). 

The characteristic low-energy feature of the calcu- 
lated spectra is a gapless mode with acoustic dispersion 
£k ~ ck for k — > 0. This remarkably reproduces a vari- 
ety of experimental results for Er 2 Ti 2 C>7 including the 
T 3 -dependence of the heat capacity [10], the inelastic 
neutron data [13], and the electron spin resonance mea- 
surements 1 1 5 j . The presence of the gapless branch in the 
harmonic approximation is directly linked to the classical 
degeneracy of the model ^ and is, therefore, accidental 
in the sense that is not related to any symmetry of the 
magnetic system. Hence, a finite gap will generically be 
produced by quantum fluctuations in higher orders of the 
spin-wave theory. Corresponding calculations are rather 
cumbersome and are not attempted here. We only re- 
mark that the generated quantum gap should depend on 
the curvature of Ai? g s . (ip) in the vicinity of its minimum. 
Since the minimum becomes more shallow with increas- 
ing anisotropy, see the inset of Fig. [2j we expect that for 
ja > 0.5 the actual gap may be rather small. 

Figure [3] shows the dispersion of four magnon branches 
in the [100] direction in comparison with the inelastic 
neutron scattering data of Ruff et al. [2] scanned from 
their Fig. 4. We find particularly good agreement be- 
tween theory and experiment for the following values of 
the microscopic parameters: J = 0.28 meV and j a = 1, 
corresponding curves are plotted by solid lines. Magnon 
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FIG. 3: (Color online) Zero-field dispersion of magnons in 
Er2Ti207 along the [100] direction. Wave-vectors are mea- 
sured in units of 2tv. Solid lines are calculated for J = 
0.28 meV and j a = 1, while dashed lines are for J = 0.33 meV 
and j a = 0.5. Open squares are experimental data points [T3] . 
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FIG. 4: (Color online) Main panel: temperature depen- 
dence of the sublattice magnetization for the classical XY py- 
rochlore antiferromagnet with in-plane anisotropic exchange 
on periodic cluster with L = 12. Inset: the fourth-order 
Binder cumulant for several lattice sizes at fixed j a = 0.5. 



energies for J = 0.33 meV and j a = 0.5 a value ex- 
tracted from the reduced moment (see above) are shown 
by dashed lines. While the low-energy branch is fitted 
equally well, a larger j a is favored by the position of the 
lower gap A « 0.38 meV. 

We can extract further information about the the mi- 
croscopic parameters from values of the critical field 
between the antiferromagnetic phase and the polarized 
paramagnetic state [T4TI17] . The high- field polarized 
states are particularly simple for field orientations along 
[100] and [110] axes. Repeating the above spin-wave cal- 
culations for H || [110] we obtain the critical field from 
the vanishing gap in the magnon spectrum 



g± f i B H[ n ^ = JS 



4 

7! 



■3a( 



73 



1 



(12) 



7.6. The dependence of (12) on 



where g± 

j a is very weak and from the extrapolated T — value 
H [ c W] ~ 1.8 T Q3H3I], we derive J ~ 0.25 meV, which is 
consistent with the estimate based on neutron data. Ad- 
ditional measurements of the excitation spectrum and/or 
the angular dependence of the critical field H c at low 
temperatures should allow a more precise identification 
of the microscopic parameters of the model (|3| in the 
case of Er 2 Ti 2 7 . 

Finite-temperature transition. — We finally consider the 
nature of ordering transition in the XY pyrochlore 
antiferromagnet with the bond-dependent exchange 
anisotropy. Since the original quantum model ([3| re- 
mains intractable for numerical investigation, we study 
instead a simplified version of ([3| with spins represented 
by classical unit vectors. The classical Monte Carlo sim- 
ulations were performed on periodic lattices of N = 4L 3 
spins with linear sizes L = 6-12 using the Metropo- 
lis algorithm with ~ 10 7 Monte Carlo steps. Temper- 
ature dependence of the sublattice magnetization m q —o 



for antiferromagnetic ordering with zero wave-vector is 
shown in Fig. [4] for three values of j a . For the model 
with isotropic exchange, j a = 0, a clear jump appears 
at T c /J — 0.110(1). However, already a small exchange 
anisotropy j a =0.1 changes the behavior m 9= o(T) to 
continuous. 

This is confirmed by the temperature dependence of 
the fourth-order Binder cumulant 



Ua = 



(13) 



According to finite-size scaling theory, Binder cumulants 
for clusters with different linear sizes L cross at the criti- 
cal point where spin-spin correlations exhibit a power-law 
decay. The crossing point T c /J — 0.445(2) for j a = 0.5 
(Fig. [4]) gives an accurate estimate for the transition point 
in this case. In our Monte Carlo simulations we also 
investigated various Potts-type order parameters, which 
distinguish between the noncoplanar ?/> 2 an d the planar 
ip3 spin structure. As T — > there is a clear selection of 
the tp2 configuration. However, at higher temperatures 
spins strongly fluctuate restoring at T — > T c the contin- 
uous XY degeneracy parameterized by angle ip in Q. 

In conclusion, we have proposed a simple quan- 
tum model for the XY pyrochlore antiferromagnet 
Er2Ti 2 07, which includes bond-dependent in-plane ex- 
change anisotropy of strength j a . We show that the order 
by disorder selection occurring for j a = is maintained 
at both quantum and classical level for perturbations up 
to at least the strength of the isotropic exchange. The ex- 
change anisotropy also reproduces the observed second- 
order transition found in Er 2 Ti 2 07 in thermodynamic 
measurements. Further work is planned to investigate 
field evolution of the antiferromagnetic state and role of 
the long-range part of the dipolar interactions, but it is 
worth noting here that the perturbation is of the same 
form, but of opposite sign to the dipolar interaction for 
n.n. pairs of spins. As noted by Stasiak et al. [53], an ex- 
change anisotropy of this kind can therefore overcome the 
n.n. part of the dipolar interaction. In addition, the long- 
range part of the dipolar interaction loses importance as 
the tetrahedral units have no net magnetic moment in 
the ground state. These observations provide a possible 
explanation as to why the dipolar effects appear to be 
absent in Er 2 Ti 2 C>7. 

The fact that our theoretical analysis with single pa- 
rameter j a is in agreement with many experimental fea- 
tures of Er 2 Ti 2 07 strongly supports the above picture 
and in the process provides strong evidence that the or- 
dering transition and the stability of the experimentally 
observed noncoplanar magnetic structure is indeed a re- 
markable experimental realization of the order by disor- 
der mechanism 12] . A consequence of this is the unique 
appearance of a non-Goldstone acoustic magnon branch 
with a very small gap. 
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Note added in proof. — After submission of this work 
we became aware of a related preprint by Savary et al. 
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